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Abstract. This paper presents the detailed results obtained in 
the search of density-velocity inhomogeneities in a volume 
limited and absolute magnitude limited sample of A-F type 
dwarfs within 125 parsecs of the Sun. A 3-D wavelet analysis 
is used to extract inhomogeneities, both in the density and 
velocity distributions. Thus, a real picture of the phase space is 
produced. Not only are some clusters and streams detected, but 
the fraction of clumped stars can be measured. By estimating 
individual stellar ages one can relate the streams and clusters 
to the state of the interstellar medium (ISM) at star formation 
time and provide a quantitative view of cluster evaporation and 
stream mixing. 

Having established this picture without assumption we come 
back to previously known observational facts regarding 
clusters and associations, moving groups or so-called super- 
clusters. 

In the 3-D position space, well known open clusters (Hyades, 
Coma Berenices and Ursa Major), associations (parts of 
the Scorpio-Centaurus association) as well as the Hyades 
evaporation track are retrieved. Three new probably loose 
clusters are identified (Bootes, Pegasus 1 and 2). The sample is 
relatively well mixed in the position space since less than 7 per 
cent of the stars belong to structures with coherent kinematics, 
most likely gravitationally bound. 

The wavelet analysis exhibits strong velocity structuring at 
typical scales of velocity dispersion a 'stream ~ 6.3, 3.8 and 
2.4 km ■ s . The majority of large scale velocity structures 
stream ~ 6.3 km ■ s" 1 ) are Eggen's superclusters (Pleiades 
SCI, Hyades SCI and Sirius SCI) with the whole Centaurus 
association. A new supercluster-liks structure is found with 
a mean velocity between the Sun and Sirius SCI velocities. 
These structures are all characterized by a large age range 
which reflects the overall sample age distribution. Moreover, 
a few old streams of ~ 2 Gyr are also extracted at this 
scale with high U components. We show that all these large 
velocity dispersion structures represent 63% of the sample. 
This percentage drops to 46% if we remove the velocity 
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background created by a smooth velocity ellipsoid in each 
structure. Smaller scales {cfstream ~ 3.8 and 2.4 km ■ s _1 ) 
reveal that superclusters are always substructured by 2 or more 
streams which generally exhibit a coherent age distribution. At 
these scales, background stars are negligible and percentages 
of stars in streams are 38% and 18% respectively. 



Key words: Techniques: image processing, Galaxy: solar 
neighbourhood - open clusters and associations - kinematics 
and dynamics - structure. 



1. Introduction 

The context of this study as well as implications of the results 
for the solar circle kinematics understanding are in Chereul et 
al (1998) (hereafter Paper II). This third paper aims to provide 
to the reader of Paper II the details of the phase space inhomo- 
geneities brought to the fore by means of wavelet analysis and 
technical details on several features of the analysis. 
Wavelet analysis is extensively used to extract deviations 
from smooth homogeneity, both in density (clustering) and 
velocity distribution (streaming). A 3-D wavelet analysis tool 
is first developed and calibrated to recognize physical inhomo- 
geneities among random fluctuations. It is applied separately 
in the position space (Section Q) and in the velocity space 
(Section ||). Once significant features are identified (whether 
in density or in velocity), feature members can be identified 
and their behaviour can be traced in the complementary space 
(velocity or density). Thus, a real picture of the phase space is 
produced. Not only are some clusters and streams detected, but 
the fraction of stars involved in dumpiness can be measured. 
Then, estimated stellar ages help connecting streams and 
clusters to the state of the ISM at star formation time and 
providing a quantitative view of cluster melting and stream 
mixing at work. 

Only once this picture has been established on a prior-less 
basis we come back to previously known observational facts 
such as clusters and associations, moving groups or so-called 
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superclusters. A sine qua non condition for this analysis to 
make sense is the completeness of data within well defined 
limits. In so far as positions, proper motions, distances, 
magnitudes and colours are concerned, the Hipparcos mission 
did care. Things are unfortunately not so simple concerning 
radial velocities and ages. More than half radial velocities are 
missing and the situation regarding Stromgren photometry 
data for age estimation is not much better. So we developed 
palliative methods which are calibrated and tested on available 
true data to circumvent the completeness failure. The palliative 
for radial velocities is based on an original combination of 



the classical convergent point method (Section 5.1) with 



the wavelet analysis. The palliative for ages (Section |3j) 
is an empirical relationship between age and an absolute- 
magnitude/colour index. It has only statistical significance in a 
very limited range of the HR diagram. 

The sampl e was pre- selected inside the Hipparcos Input 
Catalogue ( ESA, 1992 ) among the "Survey stars". The lim- 
iting magnitude is m v < 7.9 + 1.1 ■ sin | b | for spectral 



types earlier than G5 ( Turon & Crifo, 1986 ). Spectral types 
from AO to GO with luminosity classes V and VI were 
kept. Within this pre-selection the final choice was based 
on Hipparcos 'ESA, 1997) magnitude (m v < 8.0), colours 
(-0.1 < B - V < 0.6), and parallaxes (tt > 8mas). 
The sample studied (see sample named hi 25 in Paper I, 
Creze et al, 1998) is a slice in absolute magnitude of this 



selection containing 2977 A-F type dwarf stars with absolute 
magnitudes brighter than 2.5. It is complete within 125 pc 
from the Sun. 

Details concerning the wavelet analysis implementation 
through the "a trou" algorithm and the thresholding procedure 
are given in Section |2|. The palliative age determination 
method is fully described in Section |J The results of the 
density analysis in position space (clustering) are given in 
Section 0. The convergent point method is explained in Section 



5.1 and the procedures estimating spurious member s and fi eld 
stars among detected streams are given in Section 



5.2.1 



and 

5.2.2| respectively. The res ults of the velocity space analysis 
(streaming) are in Section 



5.3 



including a detailed review of 
observational facts and a comparison with evidences collected 
by previous investigators. Conclusion in Section || summarizes 
the different results. 



2. Implementation of the wavelet analysis 

Main reasons for the choice of this algorithm developed by 
Holschneider et al (1989) as well as theoretical way to compute 
wavelet coefficients are provided in Paper II. Here we focus on 
the practical computation of the wavelet coefficients (Section 
2.1) and the thresholding procedure adopted (Sections 2.2 and 
23). 
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Fig. 1. Implementation of the wavelet analysis 
2.1. The "a trou" algorithm 

The adopted mother wavelet ^{x) is defined as the difference 
at two different scales of a same scaling function $(cc) (or 
smoothing function): 



*(a0 = *(:«0 --*(-) 



(l) 



This particularity allows to use the "a trou" algorithm imple- 
mentation of the wavelet transform. This algorithm avoid the 
direct computation of the scalar product between ^(x) and the 
signal F(x) to obtain wavelet coefficients. It uses an iterative 
schema based both on the relation existing between the mother 
wavelet and the scaling function (Equation [l]) and the following 
relation 



1 



2 S + 



E 



Z=-2 



I) 



(2) 



where h(l)={jg, j, | , j, jg} is a one-dimensional discrete 
low-pass filter. h{l) is applied in the algorithm to compute iter- 
atively the different approximations (or smoothing) C s of the 
signal at each scale s. For a 3D signal like our density or veloc- 
ity distributions, h(l) is applied separately in each dimension 
to obtain the signal approximation at scale s and pixel (i,j,k): 

2 2 2 

C a (iJ,k) = E E E Kl)h(m)h(n) 

l = -2 rn=-2 n=-2 

■C s -x(i + 2 s " 1 /, j + 2 s - 1 m, k + 2 s - 1 n) (3) 

The distance between two bins increases by a factor 2 from 
scale s — 1 to scale s. The result at each scale s (C s -i(i,j, k)- 
C s (i,j, k)) is the signal difference which contains the informa- 
tion between these two scales and is the discrete set of 128 3 
wavelet coefficients W s (i,j,k) associated with the wavelet 
transform by ^(x, y, z) (Figure [l]): 

W.{i,j, k) = C._i(i,j, k) - C s (i,j, k) (4) 
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2.2. Extracting significant coefficients 

In order to remove non-zero wavelet coefficients generated by 
noise fluctuations in the 3D distributions, a local thresholding 
is applied at each scale in the space of wavelet coefficients. 
Thresholds are set at each scale and each position by estimating 
the noise level generated at the same scale by a uniform random 
signal built with the same gross-characteristics as the observed 
one at the position considered. 

- Local hard-thresholding procedure 

The noise N s (i,j, k) is estimated at the position of each 
pixel k) from the observed distribution in a filter area 
corresponding to scale s. The signal estimation S s (i,j, k), 
in this area, is: 



cients). 

The probability is 

P(W a (i,j,k)>X+) 



IwT=xj f(W s )dW s 
fw™of(W s )dW s 



(7) 



where / is the frequency of positive wavelet coefficients 
at scale s. P was set to achieve the best compromise 
between false detections and losses at each scale with a set 
of uniform background simulations containing gaussian 
structures of different sizes and signal to noise ratio. This 
optimum was found between 5-10 -4 and 10~ 4 in all 
experiments and in all scales under reasonable signal to 
noise ratio. Then e = 1CP 4 was fixed for the rest of the 
calibration. 



i+2 B j+2 a k+2 B 

S s (i,j,k)= ^2 ^2 ^2 F ( l ' m i n ) 

l=i-2 B m=j-2 B n=k-2 s 

the noise estimation is: 



(5) 



AT f ■ u\ S s(i,j,k) 
N s (i,J,k) = J (2fl+1)3 



(6) 



- Calibration curves 

A calibration curve is obtained for each scale with an 
other set of different Poisson noise simulations. Noises 
N s (i,j,k) are estimated at each pixel (i,j,k) from the 
observed distribution in the same way as described in 
The mean N. 



Section 2.2 The mean N s of Poisson noise distribution 
obtained is linked with threshold A+ (or X~) defined from 
P giving the calibration curve for scale s. 



At scale s, the number of pixels included in the convolution 
filter volume is (2 S+1 ) 3 . 

The local threshold Xf(i,j,k) is inferred from the local 
noise N s (i,j, k) through calibratio n cu rves established by 
numerical simulations (see Section 2.3) which link it to the 



noise N s (i,j,k). Then the hard-thresholding on wavelet 
coefficients follows: 



W s (i,j,k) = 



W s (i,j,k) i£W a (i,j,k)>X+(i,j,k) 







otherwise 



- Segmentation procedure 

A segmentation procedure is applied on thresholded coef- 
ficients to localize significant over-densities. This segmen- 
tation algorithm detects all maxima in a moving box of fil- 
ter size volume. Neighbouring pixels are aggregated around 
each maximum until the intensity of wavelet coefficient 
stop decreasing. This is the classical method of valley lines. 
All aggregated coefficients around a maximum are consid- 
ered to sign one structure. All structures are labelled and 
defined pixel by pixel. This procedure has the advantage of 
splitting coalescent structures. 

2.3. Thresholding calibration 

- Choice of significant probability 

Following Lega et al (1996), thresholds A are defined with 
respect to the probability P that the coefficient W s (i,j, k) 
is larger than A+ (or smaller than Xj) at a given scale s in 
the distribution of positive coefficients (or negative coeffi- 



3. Palliative age determination method 

We succeed in estimating indiv idual ages from Stromgren 
photometry ( ^igueras et al, 1991 , Masana, 1994, Asiain et al, 
1997) for a third of the sample (see Paper II). To rule out any 
bias coming from an analysis of this sub-sample alone, we pro- 
pose an empirical palliative age estimation method based on 
the Hipparcos absolute magnitude and colour for the rest of 
the sample. On a first step, we use existing ages to draw a 
plot of ages versus (Mv, (B — V)). A primary age parame- 
ter is assigned as the mean age associated to a given range of 
(Mv, (B — V)) (Section 3.1). Then Stromgren age data are 
used a second time to assign a probability distribution of pal- 
liative ages as a function of the primary age parameter (Section 



3.1. Primary age parameter 

Stromgren ages are available for 1077 A-F type stars with spec- 
tral types later than A3 (B ~ V > 0.08) because the metallic- 
ity cannot be determined for A0-A3 spectral types. With this 
sub-sample a relation between the absolute magnitude M v , 
the colour indice B — V and the age has been calibrated on 
a grid to extrapolate a less accurate age for the rest of our 
sample (1900 stars). Since the relation cannot be calibrated 
with Stromgren ages for spectral types between AO and A3, 
we add 30 known very young stars (~10 7 yr) belonging to the 
Centaurus-Crux association (see Section Q) which have B — V 
< 0.08. Consequences on the infered age distribution are in- 



vestigated in Section 3.2. Available Stromgren ages are aver- 



aged on a 200x200 grid, ranging from [-0.5,2.5] in absolute 
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Uncorrected log(AGE)^ Mv (b-v)) 



Fig. 2. Primary age parameter versus (M v , (B — V) ) relation- 
ship calibrated with the Stromgren sub-sample. Grey levels rep- 
resent the log(age) value. 



Fig. 3. Primary age parameter versus Stromgren age. 



magnitude and [-0.1,0.6] in B — V. The result is smoothed 
and extrapolated to empty cells, where feasible, using a 10x10 
moving window. This process produces a primary age param- 
eter (Figure fy. In Figure |[ Stromgren age indicators are plot- 
ted against this primary age parameter. Obviously, there is not 
a one to one relation between the primary age parameter and 
the Stromgren age. Very young Stromgren ages are the most 
affected by the degeneracy of this relation; while the bulk of 
primary age parameters older than log(age) — 8.7 are in rel- 
atively good agreement with Stromgren ages. Only 14 among 
the 1900 stars have {M v , (B - V)) out of the calibrated grid 
and do not have primary age parameter. On the basis of Figure 
[|, we can assign a probability distribution for the Stromgren 
age to each primary age parameter. 

3.2. Palliative ages 

Given a primary age parameter i, Figure || can be read as giving 
a discrete probability distribution P(j/i) of Stromgren ages j 
under i. P(j /i) is given by: 



P(j/i) = 



where n(i, j) is the number of stars in cell and 



(8) 



(9) 



These probabilities are discrete on a grid with an adaptive bin 
size to keep at least 3 stars per bin. This process produces a 
palliative age distribution which is free from possible biases 





9.5 



Fig. 4. Stromgren age distribution obtained with 1077 stars 
(full line) and palliative age distribution for the same sub- 
sample ( dotted line). 



affecting the sub-sample of stars with Stromgren photometry 
for stars with B -V > 0.08. 

The small deviations from the distribution of original 
Stromgren ages (Figure are due to finite bin steps used in 
the discretisation process. The palliative age distribution of the 
sub-sample without Stromgren photometry (Figure |]) shows 
a great difference for very young ages with respect to the 
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7 7.5 8 8.5 9 9.5 10 

log(AGE) 

Fig. 5. Stromgren age distribution obtained with 1077 stars 
(full line) and palliative age distribution for the remainder, i.e 
stars without Stromgren photometry (dotted line). 

Stromgren age distribution. The great peak at ~10 7 yr partly 
results from the poor calibration of the relation between age 
and (M V ,(B - V)) for B - V < 0.08. If the calibration of 
the relation is realized without the additional Centaurus-Crux 
stars, we would obtain a frequency, at ~10 7 yr, of ~ 170 stars 
instead of ~ 440 stars. It means that there are at least 170 
stars of few 10 7 year old and that the difference of ~ 270 stars 
is composed of A0-A3 type stars which may be older. How 
older are they? Using theoretical isochrones from Bertelli et 
al (1994), a relation between the total lifetime of a star at the 
turnoff point and the color indice (B — V) can be infered. 
We can show that the total lifetimes of A0-A3 type stars are 
between 4 - 6 10 s yr, if the metallicity of the stars are between 
Z=0.008 and Z=0.02. Taking the half lifetime of the star as the 
most probable age, we obtained, for these 270 stars, an age 
distribution in the range 2-3 10 s yr (log(age) = 8.3 — 8.5). 
This distribution of these half lifetimes remains very well 
separated from the older stars in the sample. 
To summarize, the great peak at 10 7 yr contains at least 170 
stars of ^10 7 yr and 270 stars for which ages may spread up to 
3 10 8 yr. These stars belong to the younger part of the sample. 

The palliative ages are statistical ages. Hence, they sometimes 
produce artifacts or young ghost peaks in some stream age 
distributions. Such dummy young peaks will always appear as 
the weak counterpart of a heavy peak around log(age) = 8.7. 
Nevertheless palliative ages permit to shed light on the age 
content of the phase space structures when Stromgren data are 
sparse. 



4. Clustering 



4.1. Searching for clusters 

(X,Y,Z) distributions range from -125 pc to +125 pc and are 
binned in a Sun centered orthonormal frame, X-axis towards 
the galactic center, Y-axis in the rotation direction and Z-axis 
towards the north galactic pole. The discrete wavelet analysis 
is performed on five scales: 9.7, 13.6, 21.5, 37.1 and 68.3 pc. 
These values correspond to the size of the dilated filter h at 
each scale. 

Stars belonging to each volume defined by significant coeffi- 
cients are collected. Due to the over-sampling of the signal by 
the "a trou" algorithm, some structures are detected on several 
scales. A cross-correlation has been done between all scales to 
keep the structure at the largest scale provided that there is no 
sub-structure at a lower scale (higher resolution). We control a 
posteriori the efficiency of the detection by calculating for each 
over-density the probability to find an identical concentration 
from the total volume of the sample in a uniform random 
population (see probability Pi in Table |Tj). The method also 
succeeds in finding over-densities with a high probability to 
be generated by a poissonian process. An iterative 2.5 sigma 
clipping procedure is done on tangential velocity distributions 
for each group to remove field stars and select structures with 
coherent kinematics (Table [j] and Figures |[ 0, [||] and |To| ). 



4.2. Detailed cluster characteristics 

The volumes selected by the segmentation procedure content 
10 per cent of the stars. After the 2.5 sigma clipping procedure 
on the tangential velocities, only 7 per cent are still in clusters 
or groups. Most of them are well known: Hyades, Coma 
Berenices, Ursa Major open clusters (hereafter OC1) and the 
Scorpio-Centaurus association. The following cases are the 
most interesting either because of newly discovered features 
associated to well known structures (Scorpio-Centaurus 
association) or simply because they were undetected up to now 
(Bootes and Pegasus 1 and 2). 



The Scorpio-Centaurus association shows three precisely 
limited over-densities with an elongated shape crossing 
the galactic disc (see scale 4 on Paper II, Figure 5): 
lower Cen-taurus-Crux, upper Centaurus-Lupus which 
both have identical proper motions and a third unknown 
clump at (I = 325.4°, b = -34.4°) with slightly different 
proper motions along I axis (see cluster 9 in Table [[]). 
We do not find previous mention of this extension of the 
Scorpio-Centaurus. Unfortunately very few of these stars 
have Stromgren photometry. In the case of Centaurus-Crux 
(see Figures^), among 33 stars, 3 stars having a Stromgren 
age are rather old. One of these stars is flagged as a double 
system in the Hipparcos catalogue, which could alter the 
Stromgren photometry. The two others probably do not be- 
long to the association but are not rejected by the selection 
on tangential velocities because of proper motions close to 
the association ones. However this association is known to 
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Table 1. Main characteristics of detected spatial structures after iterative 2.5sigma clipping procedure on tangential velocities: 
galactic coordinates (l,b), distance D, mean tangential velocities (Ti,Tb), tangential velocity dispersions {o~Tt,crT b ), number of 
selected stars N se ;. Other parameters are given: scales of detection S, volume v occupied by the total over-density, number of 
expected stars N e in this volume according to the mean density, number of observed stars N Q before selection on tangential 
velocities, probability Pi to find such a concentration from a poissonian distribution in this v volume and probability P2 to find 
this volume v with this concentration within the 125 pc radius sphere. A value of .001 indicates a probability P<0.001. 
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be an O-B association, which means that it is younger than 
a few tens of million years. 

The palliative age distributions show that all the three 
clumps are composed by the youngest stars of the sample. 
The analysis of the velocity field (see Section 5.3.1 ) reveals 
that these density inhomogeneities are embedded in a 
stream-like structure moving through all the 125 pc sphere. 

Two new structures are found both probably on their 
way towards disruption. In the following these loose 
clusters are designated after the constellation they are 
found in. 

2. Bootes (Figure is composed of only 5 stars at a mean dis- 
tance of 105 pc and centred on (I = 97.6°, b = +56.6°). 
This small clump is interesting for several reasons. It is 
located at the same x and z coordinates as Coma Ber. 
OC1 (Figure ^) but 62 pc ahead along y-axis with exactly 
the same space velocity components (U,V,W)=(-3.1, -7.8, 
-0.8) km • s _1 . These coincidences suggest a common ori- 
gin for these two structures. But their respective age dis- 
tributions (Figure [7] and ||) show that Bootes is probably 
younger (peak at 10 7 yr for the palliative age distribution) 



than Coma Ber. OC1 which is ~ 6 • 10 8 yr old. Moreover, 
the smaller tangential velocity dispersions of Bootes tends 
to confirm this hypothesis. A more detailed investigation 
among other spectral types is necessary to better define 
Bootes in velocity space as well as in age content. 

3. Pegasus appears to be composed of two different velocity 
components (Figure ^, |Io|). This feature do not permit 
to extract coherent velocity structures by means of the 
sigma clipping procedure. A hand selection was realized 
on the basis of the whole tangential velocity distribution. 
Pegasus 1 with 10 stars and Pegasus 2 with 8 stars are 
obtained with mean distances of respectively 97 and 89 pc 
and centred on (I = 60.7°, -32.5°) and (54.2°, -31.3°). 
Their age distributions are not very precise because of the 
lack of Stromgren ages. However, it seems that Pegasus 1 
has an age peak at 8 • 10 8 yr. It has no very young ages 
but 3 Stromgren ages are between 8 • 10 8 and 1.2 • 10 9 
yr which would mean an age older than Pegasus 1. This 
complex structure gives an instantaneous snapshot of the 
interactions affecting groups, clusters or associations: 
merging-like process can take part in dissolving spatial 
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Fig. 6. Lower Centaurus-Crux after selection on tangential ve- 
locities (top) and its age distribution (bottom). 



Fig. 7. Coma Berenices open cluster after selection on tangen- 
tial velocities (top) and its age distribution (bottom). 



inhomogeneities. 



5. Streaming 

The sub-sample with observed radial velocities is incomplete 
and contains biases since the stars are not observed at random. 
To keep the benefits of the sample completeness, a statistical 
convergent point method is developed to analyze all the stars 
(Section |j). However, this method creates spurious me mbers 
among detected streams with wavelet analysis, Section 5.2.1 
describe the procedure to handle their proportion. Moreover 
the fraction of field stars detected as stream members is also 



evaluated through the procedure in Section 5.2.2 



5.7. Recovering 3-D velocities from proper motions: a 
convergent point method 

U, V and W velocities are reconstructed from Hipparcos tan- 
gential velocities by a convergent point method for stars which 
belong to streams. All pairs of stars are considered; each pair 
gives a possible convergent point (assuming that both stars 
move exactly parallel) and the radial velocity is inferred for 
each component. Then fully reconstructed velocities of the two 
stars, V\ and V%, are considered only if | V\ — V% | does not 
exceed a fixed selection criterion. This pre-selection of recon- 
structed velocities eliminates most of false reconstructions. In 
the following, the criterion is fixed to 0.5 km ■ s . The con- 
vergent point algorithm (see Figure follows the steps: 

1. keep a pair of stars (Si, S2) 

2. define vectors A/"i and N 2 which are perpendicular respec- 
tively to the plane (OS^V^) and (OS 2 ,V 2 \\) 
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Fig. 8. Bootes loose cluster after selection on tangential veloc- Fig. 9. Pegasus 1 after selection on tangential velocities (top) 
ities (top) and its age distribution (bottom). and its age distribution (bottom). 



3. obtain JVi x = OA which is the direction of the hy- 
pothetical convergent point A 

4. test the coherence of this direction with both tangential ve- 
locities: if sgn(OA • Viy) 7^ sgn(OA ■ V2\\), it is not a 
convergent point. Go to 1 

5. calculate angles between star directions and the convergent 
point: $i=(OS7TO-4) and $ 2 =(OS7TOA) 

6. infer moduli and signs of radial velocities: 



V R = 



if $ < - 

tan $ 2 



or 



Vr = 7- — 7^ V *f $ > ? and $ ^ n 

tan($ — it) 2 



(10) 



7. calculate space velocities V\ and V2 which vectors are 
strictly parallels by construction 



V 



IV? 



(id 



8. test agreement between Vi and V2 within tolerance r\ 
I V\ — V2 \< r] with fixed 77 = 0.5 km ■ s^ 1 . 



Following this process, a large number of may-be velocities 
are calculated. Several definitions of each star velocity are ob- 
tained. All are distributed along a line in the 3D velocity space, 
part of them being spurious. Real streams produce over-density 
clumps formed by line intersections. The wavelet analysis de- 
tects such clumps in the (U, V, W) distributions. The detec- 
tion sensitivity is tested numerically by simulating a variety 
of stream amplitudes and velocity dispersions over a velocity 
ellipsoid background. Simulations show that our wavelet anal- 
ysis implementation is able to discriminate streams formed by 
at least 16 stars non-spatially localized and moving together 
with a typical velocity dispersion of 2 km ■ s^ 1 in a velocity 
background matching the sample's one. The scale at which the 
stream velocity is detected is a measure of the stream velocity 
dispersion. A more accurate knowledge of this parameter is ob- 
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Fig. 10. Pegasus 2 after selection on tangential velocities (top) 
and its age distribution (bottom). 



tain ed after the full identification of the members (see Section 

Velocities of open cluster stars are poorly reconstructed by this 
method because their members are spatially close. For such 
stars, even a small internal velocity dispersion results a poor 
determination of the convergent point. For this reason, we have 
removed stars belonging to the 6 main identified space concen- 
trations (Hyades OC1, Coma Berenices OC1, Ursa Major OC1 
and Bootes 1, Pegasus 1, Pegasus 2 groups) found in the pre- 
vious spatial analysis. Eventually, the reconstruction of the ve- 
locity field is performed with 2910 stars. 
Reconstructed (U,V,W) distributions are given in an orthonor- 
mal frame centred in the Sun velocity in the range [-50,50] 
km ■ s _1 on each component. The wavelet analysis is per- 
formed on five scales: 3.2, 5.5, 8.6, 14.9 and 27.3 km ■ s -1 . 
In the following, the analysis focuses on the first three scales 
revealing the stream-like structures because larger ones reach 
the typical size of the velocity ellipsoid. Once the segmenta- 




Fig. 11. Implementation of the convergent point method. Viu, 
V2|| are the tangential velocities measured by Hipparcos for 
two stars. V\r, V2R, are deduced radial velocities assuming 
that stars Si and S2 belong to the same stream. 



tion procedure is achieved, stars belonging to velocity clumps 
are identify. The set of velocity definitions of a star may cross 
two velocity clumps. In this case, the star is associated with the 
clump in which it appears most frequently. 

5.2. Cleaning stream statistics 

5.2.1. Estimating the fraction of spurious members in velocity 
clumps 

Despite the pre-selection of reconstructed velocities ( | Vi- V2 1 < 
rj) some spurious velocities are still present in the field. For 
this reason, real velocity clumps do include some proportion of 
spurious members generated by the convergent point method. 
Estimating the proportion of spurious members in each veloc- 
ity clump is essential and is done by comparing the mean of 
reconstructed radial velocities of each star with its observed 
radial velocity whenever this data is available (1362 among 
2910 stars). For each star in a stream, the following procedure 
is adopted: 

1. Only radial velocities reconstructed with other suspected 
members of the same stream are considered. 

2. Mean V R rec . and dispersion ay Rrec of the reconstructed 
radial velocity distribution are calculated (Figure |l2| ). 
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Fig. 12. Example ofVji rc( . distribution for one star suspected 
to be member of a stream and its observed radial velocity (dot 
dashed line ) 




Fig. 13. Distribution of weighted radial velocity deviates (nor- 
malized residual) for all suspected stream members with ob- 
served radial velocity at scale 2. Dashed line delimits the area 
of high probability of false detection. 



3. The star is confirmed to belong to the stream when the nor- 
malized residual 



R 



N 



VRobs. — V Rre 



'V Rr 



(12) 



doesn't exceed a value k. Residuals take into account er- 
rors oy Robs on the observed radial velocities given in the 
Hipparcos Input Catalogue. These errors are classified into 
4 main values: 0.5, 1.2, 2.5, 5 km - a -1 . 
The threshold k is fixed empirically on the basis of the nor- 
malized residual histogram of all suspected stream mem- 
bers with observed radial velocities (Figure |l3|) at scale 2. 
It is set to | k |=3 in order to keep the bulk of the central 
peak. The results of this paper are robust to any reasonable 
change of | K | between 2.5 and 4. 

Noise estimations in each velocity clump for the three first 
scales (provided that they have at least 3 stars with observed 
Vr), following the procedure quoted above, are shown in Fig- 
ure [l4| There are two regimes in these noise estimations: 
streams with more than 50 initial suspected members (Ni n i t ) 
which have a contamination by spurious members around 30%; 
streams with less than 50 initial suspected members which may 
have a contamination up to 85%. In this extreme case, should 
we say that these streams are false detection? Not necessar- 
ily, because our denoising method is too drastic towards small 
streams. Indeed, in those, each star has very few reconstruc- 
tions of its radial velocity with the other members of the same 
stream. The dispersion of the reconstructed Vr distribution is 
necessarily small, implying an important normalized residual. 



Then, the star is often rejected. If we refer to our previous simu- 
lations there is a high probability that under 16 initial members, 
streams are false detection. 

At the end of this selection we obtain the number of confirmed 
members among stars with observed Vr in each stream for the 
three first scales (column N se i in Tables ^, |] and |j) and their 
sum (line Tatall in column N se i). 



5.2.2. Estimating the proportion of field stars in velocity 
clumps 

The fraction of a smooth distribution filling the velocity ellip- 
soid of our complete sample, expected inside the velocity vol- 
ume spanned by the 6 superclusters described bellow, range 
between 2% and 4% depending on the position of the structure 
with respect to the distribution centroid. Adding up these con- 
tributions, 19.2% of field stars should be expected to fill the to- 
tal volume occupied by superclusters with pure random coinci- 
dence. This is about 20-30% of the stars detected as superclus- 
ter members at scale 3. However, streams with smaller velocity 
dispersions (scales 1 and 2) are not significantly affected by 
this background. Proportions of field stars for the largest struc- 
tures found at scale 3 are given in || at column % fi e id while it 
is neglected for the remaining streams. 

5.3. Stream phenomenology 

Tables § | and @ give mean velocities, velocity dispersions 
and numbers of stars remaining after correction of spurious 



members (procedure 5.2.1) and field stars (procedure 5.2.2) 
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Fig. 16. Space distribution of Hyades SCI from the Vr se- 
lected sub-sample at scale 3 (stream 3-10 in Table ^|). 



for streams at respectively scale 3, 2 and 1. Each stream has 
Ny r observed radial velocity members. Out of the Ny r stars 
with radial velocities among suspected stream members, only 
N se i get confirmed by procedure $.2.l\ So the ratio N se i/Ny r 
is an estimate of the confirmed /suspected ratio in each stream. 
Applying this ratio to Ni n u (total stream member candidates) 
we get the expected number of real stream members in each 
stream, and the total number of stream members in the sample 
(Total 3). The percentage of stars in streams in the total sample 
follows (Total 4). The correction for the uniform background 
contribution is negligible at scales 1 and 2; it is significant at 
scale 3 where the fraction of stars in streams drops from 63.0% 
to 46.4%. In the case of large velocity dispersion structures at 
scale 3 proportions of field stars is also given in column % fi e id 
and the percentage of remaining stream stars is done in column 
Ni n u line Total 5. 



velocity dispersions. 

The analysis of the age distribution inside each stream is per- 
formed on three different data sets: 

- the whole sample (ages are either Stromgren or palliative), 

- the sample restricted to stars with Stromgren photometry 
data (without selection on radial velocity), 

- the sample restricted to stars with observed (as op- 
posed to reconstructed) radial velocity data (ages are either 
Stromgren or palliative). 

The selection on photometric ages gives a more accurate 
description of the stream age content while the last sample 
permits to obtain a reliable kinematic description since stream 



members are selected through the 5.2.1 procedure. All mean 



velocities and velocity dispersions of the streams are calculated 
with the radial velocity data set. Combining results from these 
selected data sets generally brings unambiguous conclusions. 



5.3.1. Particulars on superclusters 

Streams appearing at scale 3 (cfstream ~ 6.3 km ■ s _1 ) cor- 
respond to the so-called Eggen superclusters. Four already 
known such structures are found: the Pleiades, Hyades and Sir- 
ius superclusters (hereafter SCI) and the whole Centaurus as- 
sociation. Moreover, evidence is given for one additional struc- 
ture not detected yet. The reason why this supercluster re- 
mained undetected is probably the small velocity offset with 
respect to the Sun's. At smaller scales (ff s trmm ~ 3.8 and 2.4 
km • s _1 ) superclusters split into distinct streams of smaller 



1. Pleiades SCI 

The analysis of the Pleiades SCI is realized in Paper II 
where it is found to be composed of two main streams of 
few 10 7 and 10 9 yr. 

2. Hyades SCI and NGC 1901 stream 

The velocity clump (stream 3-10 in Table ^) identified at 
scale 3 as the Hyades SCI (see Figures [TJ] for velocity and 
age distributions and Figures |l6| for space distributions) 
is located at (U,V,W)=(-32.9, -14.5, -5.6) km ■ s" 1 with 
velocity dispersions (ajj ,ay ,<rw)=(6.6, 6.8, 6.5) km ■ s _1 . 
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Fig. 15. Hyades Scl. Thresholded wavelet coefficient isocontours at W=-2.4 km ■ s _1 of the velocity field at scale 3 (top left) and 
scale 2 ("top right). This slice in wavelet coefficients at scale 2 reveals two of the three clumps composing the whole supercluster. 
Age distributions of the whole group (stream 3-10 in Table ^) at third scale (middle left) and of the first sub-stream (stream 
2-10 in Table ^|) at second scale (middle right). Age distributions of the two last sub-streams (stream 2-18 and 2-25 in Table ||) 
discovered at scale 2 (bottom). 
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The mean velocity deviates slightly from the definition 
given by Eggen (1992b) (cf. Table p|). At this resolution the 
bulk of star ages is between 4 • 10 yr and 2 • 10 9 yr with 
two peaks at 6 • 10 s and 1.6 • 10 9 yr in Stromgren age distri- 
bution plus a 10 7 yr peak in the palliative age distribution. 
Eggen (1992b) pointed out that the supercluster contains 
at least three age groups around 3 to 4, 6 and 8-10 8 yr. 
The velocity pattern splits into 3 groups at scale 2, namely 
2-10, 2-18 and 2-25 (Table ||). Each stream presents 
a characteristic age distribution, although the velocity 
separation (centers deviates from each other by several 
km ■ s^ 1 in W) does not produce a neat age separation. 
Three different main components of 10 7 , 5 — 6 • 10 8 and 
10 9 yr are mixed in the 3 clumps. The first clump peaks 
at 10 9 yr in Stromgren ages but contains a 6 • 10 s yr old 
component also revealed by palliative ages. The second 
clump peaks at 5 • 10 s yr and 10 9 yr in Stromgren ages. 
Palliative ages produce a 10 7 yr peak which is probably 
a statistical ghost of the 5 • 10 s year old c omp onent (see 
explanation of ghost at the end of Section 3.2). The third 
clump is dominated by a 5 • 10 8 yr old component with 
two older groups of 10 9 and 2 ■ 10 9 yr. One more time the 
very young peak in palliative ages is also probably due to 
the 5 ■ 10 8 year old component. 

The presence of older supercluster members around 
1.6 • 10 9 yr as stipulated by Eggen and stressed by Chen 
et al (1997) is detected in the third velocity clump. Scale 1 
does not reveal more information so that we cannot obtain 
one age for each stream. 

So, the Hyades SCI contains probably three groups of 
5 — 6 • 10 8 yr, 10 9 and 1.6 — 2 • 10 9 yr which are in 
an advanced stage of dispersion in the same velocity 
volume. Only part of these 3 streams can be linked to 
the evaporation of known open clusters. The Hyades OC1 
recent evaporation is clearly found separately in stream 
2-15. The Praesepe OC1 mean velocity (Table g) accurately 
match none of the 3 stream velocities but could explain 
the stream 2-18 despite a difference of ~ 9 km-s -1 in the 
V componen t. The NGC 1901 supercluster described in 
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Eggen, 1996| and assumed to be a Hyades SCI component 
is found separately at scale 2 (stream 2-29 in Table f| and 
exhibits a single mode in age distribution at 8 • 10 8 yr. Its 
velocity is more dissociated from the supercluster mean 
velocity than the 3 other streams which explain a best 
member extraction. 



3. SiriusSCl 

The Sirius supercluster (see Figures |l7| for velocity and 
age distributions and Figures n8[ [l^ for space distributions) 
is found on scale 3 (stream 3-19 in Table □) at mean ve- 
locity (U,V,W)=(+14.0, +1.0, -7.8) km ■ s~" with velocity 
dispersions (ct[/,0v,(7w/)=(7.3, 6.4, 5.5) km ■ s 1 . Eggen 
(1992c) identifies two age groups, 6.3 • 10 s and 10 9 yr and 
notices that there are also younger (2.5 • 10 8 yr) and older 
members (1.5 • 10 9 yr). At the coarse resolution (scale 3), 
the age distribution is in relative good agreement with this 
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Fig. 18. Space distribution of Sirius SCI from the Vr selected 
sub-sample at scale 3 (stream 3-19 in Table 



description: Stromgren ages peak at 6 • 10 8 yr and there is 
a significant proportion of stars between 10 9 and 2 • 10 9 yr. 
Stars younger than 2.5 • 10 8 yr are probably not a statistical 
ghost of the 6 • 10 8 year old component since some stars 
with Stromgren ages are also present. 
At scale 2, the supercluster splits into two distinct streams 
(stream 2-37 and 2-41 in Table ^) at respectively (U,V,W)= 
(+12.4, +0.7, -7.7) km ■ s' 1 with (cj [ /,cry,cr w )=(4.0, 4.6, 
4.7) km ■ s" 1 and (U,V,W)=(+12.4, +4.2, -9.0) km • s" 1 
with (<tu,<tv>o'w) = (3.7, 3.3, 2.9) km ■ s^ 1 producing 
a very clear age separation: the very young stars are 
separated from a part of the oldest components (6 • 10 8 
and 1.6 • 10 9 yr). The very young component appears 
exclusively in stream 2-37 (middle right of Figure [l7| ). 
At the highest resolution, on scale 1 (bottom of Figure 
|l7| ), the stream 2-41 contains oldest components still 
interpenetrated. Space distributions (Figure |l9|) show that 
the first stream, which contains the 10 7 year old component 
is still concentrated. There are too few members in the 
second clump to make conclusions. 

The Sirius SCI is composed by three age components of 
~ 10 7 , 6 • 10 8 and 1.5 • 10 9 . The younger stream is still 
concentrated both kinematically and spatially while the 
two oldest streams are mixed in a larger volume of the 
phase space. 

4. IC 2391 SCI 

The IC 2391 SCI (see Figures ^ for age distributions 
and Figure for space distributions) is not found at large 
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Fig. 17. Sirius Scl. Thresholded wavelet coefficient isocontours at W=-10.1km ■ s _1 of the velocity field at scale 3 (top left) and 
scale 2 (top right). At scale 2, Sirius SCI is composed of 2 main streams (stream 2-37 and 2-41 in Table The 2 nd is shown 
on this W slice of wavelet coefficients. Age distributions of the whole Sirius SCI at third scale (middle left) and stream 2-37 at 
second scale ( middle right). Age distributions of stream 2-41 at scale 2 (bottom left) and stream 1-56 (in Table at scale 1 
(bottom right). This latest figure (highest resolution) shows that separating oldest populations is out of reach. 
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Fig. 19. Space distributions of the 2 sub-streams of Sirius 

SCI from the Vr selected sub-sample at scale 2: stream 2-37 
(top) and 2-41 in Table [| (bottom). 



scale because it may have been merged into the Centaurus 
association velocity group. It appears separately at scale 
2 (stream 2-14 in Table g) at (U,V,W)=(-20.8, -14.5, -4.9) 
km ■ s^ 1 with velocity dispersions (ajj ,ay ,lTvk)=(4.3, 
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Fig. 20. IC 2391 SCI. Age distributions for the IC 2391 SCI 
(stream 2-14 in Table ^|) at scale 2 (top), for sub-stream 1-20 
(in Table at scale 1 (middle) and for sub-stream 1-25 (in 
Table W at scale 1 (bottom). 
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Fig. 21. Space distribution of IC 2391 SCI from the V fl se- 
lected sub-sample at scale 2 (stream 2-14 in Table [|). 



4.9, 5.0) km ■ s 1 (see Figure 15 of Paper II). Eggen 
(1991) states that IC 2391 SCI contains two ages: 8 • 10 7 
and 2.5 • 10 8 yr while Chen et al (1997) found a mean 
age of 4.6 ± 1.6 • 10 s yr. The Stromgren age distribution 
is quite different from the palliative age distribution at 
coarser resolution. Stromgren ages peak at 8 • 10 s but 
with ages up to 2 • 10 9 yr. Palliative ages exhibit a peak 
at 6 • 10 8 yr and a 10 7 year old component (Figure p0| ). 
This last peak is certainly real because its proportion is 
too high to be a statistical ghost of palliative ages from 
a 6 • 10 s year old component and moreover Stromgren 
ages show the presence of young stars. Two sub-streams 
are found at scale 1 (stream 1-20 and 1-25 in Table 
at (U,V,W)=(-20.1, -12.8, -5.0) with (au,a v ,a w )= ( 2.9, 
3.1, 1.8) km ■ s- 1 and (U,V,W)=( -20.9, -10.0, -6.1) with 
((7[/,cy,<7VK)=(4.5, 3.7, 2.9) km ■ s _1 . The stream 1-20 
contains all the youngest stars while the sub-stream 1-25 is 
only constituted of the 6 • 10 s year old population. Velocity 
dispersions of the two streams are in agreement with this 
view: they are smaller for the stream with the younger 
component. This configuration is exactly the opposite of 
the Pleiades' one: in this case the youngest population is 
more concentrated in the velocity space and is entirely 
detected in one sub-stream while the oldest span over the 
two streams. 

Centaurus Associations and the Gould belt 

Lower Centaurus-Crux and upper Centaurus-Lupus asso- 
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Fig. 22. Streams associated with Centaurus Associations. 

Age distributions of the overall association (stream 3-15 in Ta- 
ble ^) at scale 3 (top). Age distributions of Centaurus-Crux 
(stream 2-26 in Table ^) at scale 2 (middle). Age distributions 
of Centaurus-Lupus (stream 2-12 in Table ^) at scale 2 (bot- 
tom). 
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ciations (see Figures |2| for age distributions and Figures 
?? for space distributions) which are the main compo- 
nents of the entire Centaurus association are detected as 
one velocity clump at scale 3 (stream 3-15 in Table |2|) 
with (U,V,W)=(-15.2, -8.4, -8.8) km ■ s" 1 with velocity 
dispersion (<7£/,av,<7w)=(8.6, 6.7, 6.1) km ■ s _1 . Scale 3 
is a too coarse resolution and the age distributions reflect 
the overall distribution. The whole Centaurus association 
is splitted into two parts at scale 2 and does not evolve at 
scale 1 (see Paper II, Figure 15). 

At scale 2, Centaurus-Crux (stream 2-26 in Table |3|) 
and Centaurus-Lupus (stream 2-12 in Table 51) are 
identified at (U,V,W)=(-13.1, -7.9, -9.3) km • s~" with 
(au,a v ,a w )=(6.2, 6.1, 5.5) km • s" 1 and (U,V,W)=(-12.4, 
-16.5, -7.4) km ■ s~ x with (<7u,<7 V ,(j w )=(6.l, 4.6, 3.1) 
km ■ s" 1 respectively. Unfortunately, as for the Pleiades 
SCI, a lot of young stars have not Stromgren photometry. 
Stromgren age distributions peak at 6 • 10 8 yr for both 
sub-streams but palliative age distributions show the 
predominance of the very young population (10 7 yr) in 
each case. 

There is a crucial lack of radial velocities for the spatially 
clustered structures Centaurus-Crux and Centaurus-Lupus: 
one fifth of stars have observed V.r. That is why these 
space clumps are visible on space distributions when 
taking into account all the stars of the detected streams but 
disappear with the sub-sample selected on observed 
(Figure ??). At scale 2, space distributions show that stars 
of the velocity substructure identified as Centaurus-Lupus 
association belong to a disk-like structure (see XZ projec- 
tion on Figure ??) tilted with respect to the Galactic disk. 
Eigenvectors of the spatial ellipsoid are calculated. The 
two vectors associated with the largest eigenvalues allow 
to define the plane of the structure, assuming it passes 
through the Sun. The ascending node of the intersection 
between this disk-like structure and the Galactic plane 
is Iq — 317° whi ch differs slig htly from usual values 
(In ~ 275° to 290°, |Poppel, 1997f . The angle between the 
two planes is i = 18.8° in agreement with previous study 
(i = 18° - 23°) 

6. A new supercluster 

Close to the Sirius SCI in velocity space, located at the 
mean velocity (U,V,W)=(+3.6,+2.9,-6.0) km ■ s" 1 with 
velocity dispersions (er;y,ov,c7ty)=(6.8,5.0,6.3) km ■ s" 1 , 
a new massive supercluster (stream 3-18 in Table ^) is 
detected at scale 3 (see Figures ^5] for age distributions 
and Figure ?? for spatial distribution). It contains almost 
twice as many members as the Sirius SCI. None of the pre- 
viously known superclusters corresponds to this velocity 
definition. Figueras et al (1997) indicate the presence of 
a velocity structure at (U,V)=(+7,+6) which they cannot 
confirm without doubt by their analysis and interpreted it 
as a possible sub-structure of Sirius SCI with a mean age 
of 10 9 yr. We confirm the existence of a supercluster like 
structure, probably never detected before because of its 



low velocity with respect to the Sun. On a kinematics basis 
it is clearly dissociated from the Sirius SCI. 
Age distributions at coarser scale are similar to the whole 
sample ones with ages ranging from 10 7 to 2.5 • 10 9 yr. 
But at least 5 sub-streams at scale 2 (see Table ||) are 
found to form this structure. These streams show age 
distributions of relatively old components. Stream 2-27 
shows an unambiguous peak at 6 • 10 s yr with few 1.6 • 10 9 
year old stars. On the basis of velocity and age content, this 
stream could originate from the evaporation of the Coma 
OC1. Stream 2-35 has stars which are 6 • 10 8 , 10 9 and 
1.6 • 10 9 year old on the basis of Stromgren photometry but 
palliative ages exhibit only one peak at 6 • 10 s yr. Stream 
2-38 shows a peak at 5 • 10 8 yr. Stream 2-43 has Stromgren 
ages between 6 — 8 • 10 8 and few 1.6 • 10 9 year old stars 
but the palliative ages exhibit only one peak at 8 • 10 8 . 
Stream 2-44 is clearly a 10 9 year old group. All the few 
very young palliative ages in each stream are probably 
statistical ghost because very young Stromgren ages are 
never present. 

Age distributions at the highest resolution (scale 1), not 
shown here, give exactly the same results as scale 2 but the 
number of stars dramatically decreases. 
This structure has the same features as the previously 
known superclusters: a juxtaposition of several little star 
formation bursts at different epochs in adjacent cells of the 
velocity field. The correlation between velocity and age is 
not always obvious because these bursts (5 — 6 • 10 s , 8 • 10 8 
and 10 9 yr) are not very recent. As in the Hyades SCI 
case, stream velocity volumes, defined by their velocity 
dispersions, are substantially recovering. 



Implications of these results on the understanding of the super- 
cluster concept are discussed in Paper II. 



5.3.2. Outside superclusters 

- Small scale streams 

While superclusters are found to split into smaller scale 
streams most of them corresponding to well defined age, a 
number of other streams are detected only at small scales. 
Stream 2-13 in Table [| (Figure ??) is a typical example of 
such a stream. Its age distribution shows a mono-age com- 
ponent of 10 9 year old and its vertical velocity is high (W=- 
15.3 km ■ s^ 1 ). Space distribution does not fill the 125 pc 
radius sphere and Figure ?? probably shows on (X,Z) and 
(Y,Z) projections the orbit tube in which stars are confined. 
Stream 2-7 in Table || shows also the same features with 
a 6 • 10 8 year old component and a lower vertical velocity 
(W=-8.5 km ■ s- 1 ). 

- Particulars on oldest groups 

Another striking feature of the velocity field is the 
existence of a 2 • 10 9 year old population still in velocity 
structures. These streams are only detected at the coarser 
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Fig. 25. New supercluster. Age distributions of the new moving group (stream 3-18 in Table ^|) at scale 3 (top left) and the 
sub-stream 2-27 (in Table ^|) at scale 2 (top right). Age distributions of the sub-streams 2-35 (middle left) and 2-38 (middle 
right ) at scale 2. Age distributions of the sub-streams 2-43 (bottom left) and 2-44 (bottom right) at scale 2. 
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resolution (scale 3) in agreement with an intrinsically 
large velocity dispersion. They have much less members 
(between 20 and 30 initial members) than the previously 
investigated streams and have very few observed Vr. 
Two main old groups (stream 3-9 and 3-14 in Table [| 
and Figures ??, ?? ) are clearly detected and have similar 
characteristics: 

- Age distributions peak between 1.6 and 2 • 10 9 yr 

- Their velocity dispersion obtained from the few se- 
lected stars on observed Vr are of order 6 km • s _1 
(stream 3-14 seems to have a lower velocity dispersion 
but the result is obtained for only 3 stars). 

- U-component is positive (towards the galactic center) 
and often larger than 20 km • s _1 

- Space distributions may still be clumpy 

We should notice at this point that these old star streams 
are both affected by two biases with respect to the whole 
sample characteristics. They have much less observed 
radial velocities (16% vs. 46%) and have much more 
Stromgren photometry (57% vs. 37%). For this reason, 
space distributions are displayed without radial velocity 
selection. 

Streams are rather fragile structures originating from 
cluster evaporation. The survival of these relatively old 
moving groups with coherent ages can be explain by 
two non exclusive mechanisms. On one hand, recent 
simulations (Zwart et al, 1998) of heavy cluster dynamical 
evolution (with ~ 32000 particles) have shown that 
contrary to lighter open clusters whose typical lifetime is 
1 - 2 • 10 8 (|Wielen, 197l| and[Lynga, 1982|), heavy clusters 



can survive up to 4 Gyr in the Galactic gravitational field. 
This long lifetime can authorize to maintain old streams 
but such heavy open clusters are probably very rare in the 
disc. On the other hand, as sketched in Dehnen (1998) to 
explain their eccentric orbits, stars of these moving groups, 
provided that they were formed in the inner part of the 
disc, could have been trapped into resonant orbits with 
the non-axisymmetric force created by the gravitational 
potential of the Galactic bar. This latest explanation does 
not require any minimum lifetime to the initially bound 
structure from which streams originate. 



6. Conclusions 

A systematic multi-scale analysis of both the space and 
velocity distributions of a thin disc young star sample has been 
performed. Over-densities in space distributions are mainly 
well known open clusters (Hyades OC1, Coma OC1 and Ursa 
Major OC1) and associations of stars (Centaurus-Crux and 
Centaurus-Lupus). Evidence is given for three nearby loose 
clusters or associations not detected yet (Bootes and Pegasus 1 
and 2). Evaporation of relatively massive stars (1.8 M©) out of 
the Hyades open cluster is also mapped in Paper II and form 



an asymmetric pattern. The sample is well mixed in position 
space since no more than 7% of the stars are in concentrated 
clumps with coherent tangential velocities. 
The 3D velocity field reconstructed from a statistical conver- 
gent point method exhibits a strong structuring at typical scales 
of ~5 stream ~ 6.3, 3.8 and 2.4 km ■ s . At large scale (scale 3) 
the majority of structures are Eggen's superclusters (Pleiades 
SCI, Hyades SCI and Sirius SCI) with the whole Centaurus 
association embedded in a larger kinematic structure. A new 
supercluster-like structure is found with a mean velocity 
between the Sun and Sirius SCI velocities. These large scale 
velocity structures are all characterized by a large age range 
which reflects the overall sample age distribution. Moreover, 
few old streams of ~ 2 Gyr are also extracted at this scale with 
high U components towards the Galactic center. Taking into 
account the fraction of spurious members, evaluated with an 
observed radial velocity data set, into all these large velocity 
dispersion structures we show that they represent 63% of 
the sample. This percentage drops to 46% if we remove the 
velocity background created by a smooth velocity ellipsoid 
in each structure. Smaller scales (Wstream ~ 3.8 and 2.4 
km ■ s _1 ) reveal that superclusters are always substructured 
by 2 or more streams which generally exhibits a coherent age 
distribution. The older the stream is the more difficult the age 
segregation between close velocity clumps is. At scale 2 and 
1, background stars are negligible and percentages of stars in 
streams, after evaluating the fraction of spurious members, are 
38% and 18% respectively. 
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Table 2. Main characteristics of detected streams at scale 3 with filter size of 8.6 km ■ s _1 (a 'stream ~ 6.3 km ■ s^ 1 ). U ± <r u , 
V ± cr v and W ± <j w are the mean velocity components and their dispersions calculated, after selection on Vr, with the true 
radial velocity data - Ni„u is the initial number of stars belonging to the structure - Ny R is the number of stars with observed 
radial velocity - N se ; is the number of stars among Ny r after selection 5.2.1. The selection has not been done in case where 
there is less than 3 observed radial velocities. In this last case and also when none of the stars are selected we do not present the 
line in the table. However the same stream identifiers are conserved: the first digit is the scale followed by the stream number. 
Column % fi e id is the percentage of field stars estimated in each structure with the procedure 5.2.2. Cross-identification is done 
with Eggen 's superclusters ( SCI) and open cluster ( OCl) data from Palous and Eggen ( in Gomez et al, 1990) ( see Table ^|) - Total 
1 gives the sum of suspected stream members (col. Ni n u) and the sum of observed Vr among them (col. Nv R ) - Total2 gives 
the sum of all confirmed stream stars among available observed radial velocities (col. N se l) - Total 3 is the sum of all expected 
stream members in each stream ( col. Ni n a ) taking into account the confirmed/suspected ratio obtained in each stream - Total 4 
gives the inferred fraction of stars in stream in the full sample of 2910 stars (col. Ni n u). Total 5 is the fraction of stars in stream 
corrected for field contamination (col. Ni n u). 
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Table 3. Main characteristics of detected velocity structures at scale 2 with filter size of 5.5 km ■ s~ 1 (a 'stream ~ 3.8 km ■ s^ 1 )- 
Legend is the same as Table ^| 
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2-40 


16 


8 


3 


-4.8± 3.1 


5.2± 6.7 


-11. 3± 3.6 


2-41 Sirius SCI (2) 


42 


21 


13 


12.4± 3.7 


4.2± 3.3 


-9.0± 2.9 


2-43 New SCI (4) 


77 


46 


28 


5.3± 2.7 


6.5± 3.6 


-11.4± 4.0 


2-44 New SCI (5) 


45 


20 


8 


6.3± 2.8 


6.5± 2.6 


0.8± 2.7 


2-46 


24 


13 


7 


-4.8± 6.1 


10.8± 2.0 


-6.5± 1.9 


Total 1 


1706 


827 










Total2 






538 


(Vr confirmed stream members) 


Total3 


1114.3 






(Total stream members) 





Total4 



38.3% 



(Fraction of stars in streams) 
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Table 4. Main characteristics of detected velocity structures at scale 1 with typical size of 3.2 km ■ s _1 ("a stream ~ 2.4 km ■ s~ x ) 
- Legend is the same as Table ^| 



Stream 


vr 


NT 




U ±a u 


T T 1 _ 

V ±a v 


W ±<7 W 








(I « 1=3.) 


(Km ■ s ) 


(km ■ s ) 


(km ■ s ) 


1-1 


14 


1 


3 


-14.9± 0.5 


-27. 8± 0.6 


-1.7± 1.9 


1-2 Pleiades OC1 


40 


17 


12 


-4.7± 2.0 


-23. 8± 1.2 


-8.9± 2.6 


1-4 


16 


6 


3 


-15. 1± 3.1 


-24. 1± 2.9 


-0.7± 1.7 


1-6 Pleiades SCI (1) 


33 


14 


9 


-13. 1± 3.1 


-21.9± 3.3 


-7.1± 2.5 


1-7 Pleiades SCI (2) 


46 


30 


20 


-11. 1± 1.7 


-21.9± 3.0 


-6.0± 1.8 


1-10 


22 


10 


3 


-31.0±5.2 


-19.0± 1.5 


-6.0± 2.7 


1-11 


18 


6 


3 


-19.2± 3.4 


-19.0± 2.4 


-7.1± 3.5 


1-12 


22 


11 


3 


-9.9± 0.4 


-20.3± 1.6 


-0.5± 1.8 


1-13 


22 


11 


3 


-16.1± 2.0 


-20.9± 0.1 


1.9± 2.9 


1-15 


32 


17 


9 


-18. 8± 4.2 


-15.7± 3.4 


-1.8± 2.3 


1-17 Cen-Lup 


61 


17 


5 


-9.4± 3.3 


-18.2± 2.7 


-5.9± 1.1 


l-19Hyades SCI (1) 


24 


12 


10 


-29.0± 5.5 


-12.5± 6.0 


-0.0± 4.4 


1-20 IC 2391 SC1(1) 


33 


15 


8 


-20.1± 2.9 


-12.8± 3.1 


-5.0± 1.8 


1-21 


8 


6 


3 


-19.9± 1.1 


-14.8± 3.2 


-15.5± 0.5 


1-22 


28 


12 


3 


-27.2± 2.4 


-17.1± 2.8 


-8.4± 3.2 


1-23 Hyades SCI (2) 


39 


22 


8 


-33.2± 2.3 


-13.1± 3.5 


-6.6± 2.5 


1-24 


32 


12 


3 


-13. 8± 2.7 


-15.5± 5.9 


1.1± 3.4 


1-25 IC 2391 SCI (2) 


41 


21 


8 


-21.0± 4.5 


-10.0± 3.7 


-6.1±2.9 


1-26 


23 


11 


6 


-2.4± 1.6 


-11.7± 1.8 


-4.5± 3.2 


1-27 


12 


7 


3 


-30.5± 1.0 


-10.8± 0.6 


-9.1± 2.1 


1-29 Cen-Crux 


45 


17 


6 


-13. 1± 3.2 


-9.4± 3.7 


-7.7± 4.5 


1-31 


53 


21 


12 


-19.5±3.4 


-5.8± 2.3 


-10.5± 3.0 


1-32 


17 


7 


4 


-25. 8± 2.0 


-9.9± 3.2 


-0.2± 3.0 


1-34 


24 


11 


7 


-8.0± 1.0 


-3.8± 1.3 


-13.2± 2.2 


1-37 


14 


9 


3 


3.1± 5.4 


-1.1± 1.8 


-15.2± 3.4 


1-39 


13 


6 


5 


22.1± 1.5 


-1.1± 1.9 


-0.5± 2.2 


1-41 New SCI (2) 


39 


18 


11 


4.1± 4.2 


0.0± 2.4 


-7.7± 1.8 


1-45 


14 


5 


3 


8.5± 1.4 


0.9± 1.3 


-2.4± 2.5 


1-46 


27 


13 


4 


1.4± 4.3 


1.4± 2.0 


-3.4± 2.0 


1 AQ 


1 u 


f, 
u 


-i 
j 


1 J .UZL 1 .D 




1 4- 1 fi 


1-51 Sirius SCI (2) 


16 


7 


4 


8.9± 3.1 


4.8± 1.6 


-12.9± 0.6 


1-54 


18 


11 


5 


6.5± 1.2 


5.9± 1.8 


-10.6± 2.1 


1-55 New SCI (3) 


18 


9 


5 


2.9± 4.4 


3.9± 1.8 


-4.9± 2.6 


1-56 Sirius SCI (1) 


13 


8 


4 


12.0± 2.7 


5.9± 1.3 


-9.2± 1.1 


1-58 


20 


11 


6 


3.7± 1.9 


7.2± 3.5 


-10.0± 2.7 


1-59 New SCI (5) 


24 


9 


4 


6.2± 0.6 


6.0± 1.2 


-0.3± 1.9 


1-61 New SCI (4) 


15 


5 


4 


7.2± 1.8 


9.9± 2.2 


-10.5± 2.7 


1-62 


17 


9 


5 


-3.9± 7.6 


10.5± 2.2 


-7.2± 1.1 


1-59 New SCI (5) 


24 


9 


4 


6.2± 0.6 


6.0± 1.2 


-0.3± 1.9 


1-61 New SCI (4) 


15 


5 


4 


7.2± 1.8 


9.9± 2.2 


-10.5± 2.7 


1-62 


17 


9 


5 


-3.9± 7.6 


10.5± 2.2 


-7.2± 1.1 


Total 1 


1262 


556 










Total2 






220 


(Vr confirmed stream members) 


Total3 


521.2 






(Total stream members) 





Total4 



17.9% 



(Fraction of stars in streams) 
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Table 5. Cross-identification data for known kinematical groups: Open clusters (OCl) and Superclusters (SCI) 



U V W Age Source 

(km ■ s _1 ) (yr) 



OCl 



Hyades 


-44.4 


-17.0 


-5.0 


8±2-10 8 


Pallousetal, 1977 


Pleiades 


-5.8 


-24.0 


-12.4 


5 ±M0 7 


Pallous et al, 1977 


IC2391 


-18.3 


-13.5 


-5.9 


3 ±M0 7 


Pallousetal, 1977 


U.Ma 


+14.5 


+2.5 


-8.5 


3 -10 8 


Eggen, 1973) 


Coma Ber. 


-1.8 


-8.2 


-0.7 


4±M0 8 


Pallous et al, 1977 


IC2602 


-0.7 


-25.7 


-1.4 


1 ±0.5-10 7 


Pallous et al, 1977 


Praesepe 


-37.1 


-23.5 


-7.0 


7 ±2-10 8 


Pallous et al, 1977 


Alpha Per. 


-10.8 


-20.5 


-0.7 


4±0.5-10 7 


Pallous et al, 1977 



SCI 



Hyades 
Hyades 


-40.0 
-44.4 


-17.0 
-17.8 


-2 .0 
-1.5 


Pleiades 
Pleiades 


-11.6 
-10.0 


-20.7 
-19.0 


-10.4 
-8.1 


IC 2391 
IC 2391 
IC 2391 


-22.4 
-20.8 
-15.9 


-17.5 
-15.9 
-13.1 


-9.4 
-8.3 
-4.5 


Sirius 
Sirius 


+15.0 
+13.0 


+1.0 

+3.2 


-11.0 
-7.5 


NGC 1901 


-26.4 


-10.4 


-1.5 



3- 6- 8-10* 
13.3-10 8 



0.7-10 8 



2.8-10 8 



11.0 6.3-10* &10 H 



6.5-10 



.8 



Eggen, 1992b 



Chen et al, 1997 



Eggen, 1992. 
Chen et al, 1997 



Eggen, 1991| 
Eggen, 1992d| 
Chen et al, 1997 



Eggen, 1992c 
Chen et al, 1997 



1.5 8-10 8 & 1.2-10 9 Eggen, 199€ 
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